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We study a family of high order Ehrlich-type methods for approximating all zeros 
of a polynomial simultaneously. Let us denote by the famous Ehrlich method 
(1967). Starting from Kjurkchiev and Andreev (1987) have introduced 
recursively a sequence (7’^^^)“=i of iterative methods for simultaneous finding 
polynomial zeros. For given > 1, the Ehrlich-type method has the order 

of convergence 2N -F 1. In this paper, we establish two new local convergence 
theorems as well as a semilocal convergence theorem (under computationally 
verifiable initial conditions and with a posteriori error estimate) for the 
Ehrlich-type methods Our first local convergence theorem generalizes a 

result of Proinov (2015) and improves the result of Kjurkchiev and Andreev 
(1987). The second local convergence theorem generalizes another recent result 
of Proinov (2015), but only in the case of maximum-norm. Our semilocal 
convergence theorem is the first result in this direction. 
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1 Introduction 

Throughout this paper (K, | • |) denotes an algebraically closed field and K[z] denotes 
the ring of polynomials (in one variable) over IK. For a given vector x in K", Xi always 
denotes the fth coordinate of a;. In particular, if F is a map with values in K", then 
Fi{x) denotes the fth coordinate of the vector F{x). We endow the vector space K" 
with a norm ||a:||p defined as usual: 


ll^llp = 



(1 < P < oo) 


and with coordinate-wise ordering ^ defined by 


X F y if and only if Xi < yi for each f = 1,..., n (1-1) 

for x,y G M". Then (M”, || • ||p) is a solid vector space. Also, we endow IK" with the 
cone norm || • ||: K" —)■ K" defined by 


||a;|| = (|a;i|,...,|a;„|). 
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Then (K”, || • ||, is a cone nornied space over R” (see, e.g., Proinov [1]). 

Let / € IK[z] be a polynomial of degree n > 2. A vector ^ G K" is said to be a 
root-vector of / if f{z) = Oq ~ Ci) for all 2 : G K, where Oq € K. We denote 

with sep(/) the separation number of / which is defined to be the minimum distance 
between two distinct zeros of /. 

1.1 The Weierstrass method and Weierstrass correction 

In the literature, there are a lot of iterative methods for finding all zeros of / 
simultaneously (see, e.g., the monographs of Sendov, Andreev and Kjurkchiev [2], 
Kjurkchiev [3], McNamee [4] and Petkovic [5] and references given therein). In 1891, 
Weierstrass [6] published his famous iterative method for simultaneous computation 
of all zeros of /. The Weierstrass method is defined by the following iteration 


^.(/c+i) ^ ^(fe) _ fc = 0,1, 2,..., 


( 1 . 2 ) 


where the operator Wf. D C K” —)■ K” is defined by 



where oq G K is the leading coefficient of / and the domain 2) of W is the set of all 


vectors in K" with distinct components. The Weierstrass method (1.2) has second- 
order of convergence provided that all zeros of / are simple. The operator Wf is 
called Weierstrass correction. We should note that Wf plays an important role in 
many semilocal convergence theorems for simultaneous methods. 

1.2 The Ehrlich method 

Another famous iterative method for finding simultaneously all zeros of a polynomial 
/ was introduced by Ehrlich [7] in 1967. The Ehrlich method is defined by the 
following fixed point iteration: 


a;('=+i) =T(x('=)), fc = 0,1,2,..., 


(1.4) 


where the operator T: S’ C K" —)■ K" is defined by T{x) = {Ti{x), ..., Tn{x)) with 


Ti{x) = Xi- 






(1.5) 



1 


and the domain of T is the set 



( 1 . 6 ) 


Here and throughout the paper, we denote by the set of indices 1,..., n, that is 
In = {1, • ■ • The Ehrlich method has third-order of convergence if all zeros of 
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/ are simple. The Ehrlich method was rediscovered by Abert [8] in 1973. In 1970, 
Borsch-Supan [9] introduced another third-order method for numerical computation 
of all zeros of a polynomial simultaneously. In 1982, Werner [10] has proved that 
the both methods are identical. The Ehrlich method (1.4) is known also as “Ehrlich- 
Abert method”, “Borsch-Supan method” and “Abert method”. 

Recently, Proinov [11] obtained two local convergence theorems for Ehrlich 
method under different types of initial conditions. The first one generalizes and 
improves the results of Kyurkchiev and Tashev [12, 13] and Wang and Zhao [14, 
Theorem 2.1]. The second one generalizes and improves the results of Wang and 
Zhao [14, Theorem 2.2] and Tilli [15, Theorem 3.3]. 

Before we state the two results of [11], we need some notations which will be used 
throughout the paper. For given vectors x G K" and y G K”, we define in K" the 
vector 


_ /kl| \^n\ \ 

y V J/i ’ ' ’ 2/" / ’ 

provided that y has no zero components. Given p such that 1 < p < oo, we always 
denote by q the conjugate exponent of p, i.e. q is defined by means of 

1 < 9 < oo and l/p-|-l/g = l. 

In the sequel, we use the function d: K" —>• K" defined by d{x) = (di(x), ..., dn(x)) 
with 


di{x) = mh\\xi — Xj\ (f = l,...,n). 


Let o > 0 and 6 > 1. We define the real function (j) by 


^ (1 — t){l — bt) — at^ 

and the real number R as follows 


5 -l-1 -l- \/(5 — 1)^ -l- 8a 


(1.7) 


( 1 . 8 ) 


Theorem 1.1 (Proinov [11]) Let f G K[z] be a polynomial of degree n >2 which 
has only simple zeros, f G K" be a root-vector of f and 1 <p < oo. Suppose 
x^^'^ G K" is an initial guess satisfying 


E{x^°^) 


^ 

diO 


<R = 


p 


2 

r+l +v^l)2+8a’ 


(1.9) 


where a = {n — 1)^/'^ and b = 2fl'^. Then Ehrlich iteration (1.4) is well-defined and 
converges cubically to f with error estimates 


||^(fe+i) _^|| ^ A3''||^(fc) _^|| iia-W _^|| ^ y(3''-i)/2||^(o) _^|| 

for all k> 0, where A = (f>{E{x^^'^)) and the function (j) is defined by (1.7). 
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Theorem 1.2 (Proinov [11]) Let f G K[ 2 :] be a polynomial of degree n> 2, f G K" 
be a root-vector of f and 1 < p < oo. Suppose G K" is a vector with distinct 
components satisfying 


E{x^°^) 


f 

d{x^°l) 


<R = 

P 


2 

6+1 +v^l)2+8a’ 


( 1 . 10 ) 


where a = {n— 1)^/'^ and b = 2^/'^. Then f has only simple zeros in IK and Ehrlich 
iteration (1.4) is well-defined and converges to f, with error estimates 

||^(fe+i) _ ^11 ^ - Cll and \\x^'^^ - ^|| ^ - ^|| 


for all k > 0, where A = (j){E{x^^^)), 9 = tp{E{x^^'^)) and the function (j) is defined 
by (1.7) and the function if 


ifit) 


(1 — 6)(1 — bt) — af^ 
1 — t — at'^ 


( 1 . 11 ) 


Moreover, the method converges cubically to f provided that E{x^^'^) < R. 


1.3 A family of high-order Ehrlich-type methods 

In the following definition, we define a sequence (T*'^^)^^o iteration functions in 
the vector space K". In what follows, we define the binary relation ff on K" by 


xffv ^ for all i,j G In with i j. (1-12) 

Definition 1.3 Let / G IK[z] be a polynomial of degree n > 2. Define the sequence 
(T(^))^^g of functions : Ejy C K” — >■ K” recursively by setting T^^\x) = x 
and 


rp{N+l) 


(x) = Xi- 


_ fjXi) 

/'(x,) - /(x*)^ — 



(f = l,...,n). 


(1.13) 


where the sequence of domains (Djv)“=o defined recursively by setting 

Dq = IK" and 


Dn-\-i = {x G Dn : xff (x), f'{xi) - /(xi) ^0 for f S /„}. 

’{x) 


(1.14) 


Given G N, the A^th method of Kjurkchiev-Andreev’s family can be defined by 
the following fixed-point iteration: 

xl^’+i) =tW(x('’’)), /c = 0,1,2,.... (1.15) 

It is easy to see that in the case A^ = 1 the Ehrlich-type method (1.15) coincides 
with the classical Ehrlich method (1.4). The order of convergence of the Ehrlich-type 
method (1.15) is 2N + 1. 
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Kjurkchiev and Andreev [16] established the following convergence result for the 
Ehrlich-type methods (1.15). This result and its proof can also be found in the 
monographs of Sendov, Andreev and Kjurkchiev [2, Section 19] and Kjurkchiev ]3, 
Chapter 9.2]). 


Theorem 1.4 (Kjurkchiev and Andreev [16]) Let f G C]z] be a polynomial of 
degree n>2 which has only simple zeros, ^ G C" be a root-vector of f and N >1. 
Let 0 < ft. < 1 and c > 0 be such that 


6 > 2c(l -I- (2n 


l)ft) and 


9 

nc 

{6 — c)(S — 2c— 2ch) — 3(n 


I)c2ft2 - 


(1.16) 


where 6 = sep(/). Suppose G C" is an initial guess satisfying the condition 


- ^lloo < eft. 


(1.17) 


Then the Ehrlich-type method (1.15) converges to f with error estimate 


— ^jjoo < for all k >0. 


(1.18) 


1.4 The purpose of the paper 

In this paper, we present two new local convergence theorems as well as a semilocal 
convergence theorem (under computationally verifiable initial conditions and with 
a posteriori error estimate) for Ehrlich-type methods (1.15). Our first local conver¬ 
gence result (Theorem 4.6) generalizes Theorem 1.1 (Proinov [11]) and improves 
Theorem 1.4 (Kjurkchiev and Andreev [16]). Our second local convergence result 
(Theorem 5.4) generalizes Theorem 1.2 (Proinov [11]), but only in the case p = oo. 
Furthermore, several numerical examples are provided to show some practical ap¬ 
plications of our semilocal convergence result. 

2 A general convergence theorem 

Recently, Proinov [17, 18, 19] has developed a general convergence theory for iter¬ 
ative processes of the type 


Xk+i = Txk, fc = 0,1,2,..., 


( 2 . 1 ) 


where T: D C X ^ X is axi iteration function in a cone metric space X. In order to 
make this paper self-contained, we briefly review some basic definitions and results 
from this theory. 

Throughout this paper J denotes an interval on K_|. containing 0. For an integer 
fc > 1, we denote by Sk{t) the following polynomial: 

Sk{t) = l + t + ...+t^-\ 

If fc = 0 we assume that Skft) = 0. Throughout the paper we assume by definition 
that 0° = 1. 



Proinov and Vasileva 


Page 6 of 26 


Definition 2.1 ([18]) A function (p: J ^ IR+ is called quasi-homogeneous of de¬ 
gree r > 0 on J if it satisfies the following condition: 

(p(t) for all A S [0,1] and t € J. (2.2) 

If m functions ipi,..., ipm are quasi-homogeneous on J of degree ri,... ,rm, then 
their product is a quasi-homogeneous function of degree ri -I- ... -I- on 

J. Note also that a function (p is quasi-homogeneous of degree 0 on J if and only it 
is nondecreasing on J. 

Definition 2.2 ([17]) A function p: J ^ J is said to be a gauge function of order 
r > 1 on J if it satisfies the following conditions: 

(i) ip is quasi-homogeneous of degree r on J; 

(ii) p{t) <t for all t G J. 

A gauge function p of order r on J is said to be a strict gauge function if the 
inequality in (ii) holds strictly whenever t S J\{0}. 

The following is a sufficient condition for a gauge function of order r. 

Lemma 2.3 ([18]) Ifp'- J —t K+ is a quasi-homogeneous function of degree r > 1 
on an interval J and R > 0 is a fixed point of p in J, then p is a gauge function 
of order r on [0, R]. Moreover, if r > 1, then function p is a strict gauge of order r 
on J = [0, R). 


Definition 2.4 ([17]) Let T: D C X —>■ X he a map on an arbitrary set X. A 
function E: D ^ M+ is said to be a function of initial conditions of T (with a 
gauge function p on J) if there exist a function p: J ^ J such that 

E(Tx) < p{E(x)) for all a: € D with Tx G D and E{x) G J. (2-3) 

Definition 2.5 ([17]) Let T ■. D <G X ^ X he a map on an arbitrary set X, and 
let E: D ^ K+ be a function of initial conditions of T with a gauge function on 
J. Then a point x G D is said to be an initial point of T (with respect to E) if 
E{x) G J and all of the iterates T^x {k = 0,1, 2,...) are well-defined and belong to 

D. 

The following is a simple sufficient condition for initial points. 

Theorem 2.6 ([18]) Let T: D C X ^ X be a map on an arbitrary set X and 
E: D ^ K+ be a function of initial conditions of T with a gauge function p on J. 
Suppose that x G D with E{x) G J implies Tx G D. Then every point Xq G D such 
that E{xo) G J is an initial points of T. 

Definition 2.7 ([19]) Let T: D C X ^ X he an operator in a cone normed space 
{X, II • II) over a solid vector space {Y, ^), and let A: D —^ IR+ be a function of initial 
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conditions of T with a gauge function on an interval J. Then the operator T is said 
to be an iterated contraction at a point ^ € D (with respect to E) if E(^) € J and 

\\Tx — ^11 ^ j3{E{x))\\x — ^11 for all x G D with E{x) S J, (2.4) 

where the control function /3: J —>■ [0,1) is nondecreasing. 

The following fixed point theorem plays an important role in our paper. 

Theorem 2.8 (Proinov [19]) Let T: D C X ^ X be an operator of a cone normed 
space {X, II • II) over a solid vector space and let E: D ^ K_|_ be a function 

of initial conditions of T with a gauge function ip of order r > 1 on an interval 
J. Suppose T is an iterated contraction at a point f, with respect to E with control 
function fd such that 

t /3(t) is a strict gauge function of order r on J (2-5) 

and there exist a function "0: J —>■ K+ such that 

j3{t) = (j)(f)il){t) for all t G J, (2-6) 

where ^: J —>■ M+ is a nondecreasing function satisfying 

ip(t) = t(j){t) for all t G J. (2-7) 

Then the following statements hold true. 

(i) The point is a unique fixed point of T in the set U = {x G D : E{x) G J}. 

(ii) Starting from each initial point x^^'^ of T, Picard iteration (2.1) remains in 
the set U and converges to ^ with error estimates 

||^(fe+i)_^|| ^ ||x('=)-C|| and ||x('=)-e|| ^0'=A^'=W||x(°)-e|| (2.8) 

for all k> 0, where A = (f){E{x^^'^)) and 9 = ’tp{E{x^^'^)). 

In the case (3 = (j), Theorem 2.8 reduces to the following result. 

Corollary 2.9 ([19]) Let T: D G X ^ X be an operator in a cone normed space 
(X, II • II) over a solid vector space (T, ^), and let E: H —> M+ be a function of 
initial conditions of T with a strict gauge function ip of order r > 1 on an interval 
J. Suppose that T is an iterated contraction at a point ^ with respect to E and with 
control function (j) satisfying (2.7). Then the following statements hold true. 

(i) The point ^ is a unique fixed point of T in the set U = {x G D : E{x) G J}. 

(ii) Starting from each initial point x^^'^ ofT, Picard iteration (2.1) remains in U 
and converges to ^ with order r and error estimates 

||^(fc+i)_^|| ^ A"" ||x('=)-C|| and ||x('=)-C|| ^ ||x(°)-e||, (2.9) 


for all k> 0, where A = 4>{E{x^^'>)). 
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3 Some inequalities in 

In this section, we present some useful inequalities in K" which play an important 
role in the paper. 


Lemma 3.1 ([20]) Let u,u€lK", v be a vector with distinct components and 
1 < p < oo. Then for all i,j G In, 


I Mi 


'^^■1 > 


1 - 


u — V 
d{v) 





(3.1) 


^J\ — 


( 1 - 2i/« 

u — V 


d{v) 


Vi - Vi 


(3.2) 


Lemma 3.2 ([19]) Let u,v G K" and 1 < p < oo. // the vector v has distinct com¬ 
ponents and 


u — V 
d{v) 


< 


p 


1 

2 


then the vector u also has distinct components. 


Lemma 3.3 ([21]) Let u,u,^gK", he a vector with distinct components, 
1 < p < oo and 

iiM-a^iiM-eii. (3.3) 


Then for all i^j G In} 




f 1 - 2i/« 


1 

d{li} 


IC. - 




(3.4) 


Lemma 3.4 Let u,v,^ G K", a > 0 and l<p<oo. If v is a vector with distinct 
components such that 


||m-CII ^ a||M-CII, 


(3.5) 


then for all i,j G In, 


-v^\> 1 - (1 + a) 




div) 


^J\' 


(3.6) 


Proof By the triangle inequality of cone norm in K" and (3.5), we obtain 


|| m - m || ^ \\u-^\\ + || m -^|| ^ (l + a)||z;-C||, 
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which yields 


u — V 
d{v) 


^ (1 + a) 


d{v) 


Taking the j9-norm, we get 


u — V 
d{v) 


<{l + a) 

P 


d{v) 


(3.7) 


From (3.1) and (3.7), we obtain (3.6) which completes the proof. 


□ 


Lemma 3.5 Let u,v,^ € K", a > 0 and l<p<oo. If v is a vector with distinct 
components such that (3.5) holds, then for all i,j S In, 


-Uj\ > l-2^/«(l + a) 


v-f 


d{v) 


(3.8) 


Proof From (3.2) and (3.7), we get (3.8) which completes the proof. 


□ 


4 Local convergence theorem of the first type 

Let / € K[z] be a polynomial of degree n> 2 which has only simple zeros in K, 
and let ^ S K" be a root-vector of /. In this section we study the convergence of 
the Ehrlich-type methods (1.15) with respect to the function of initial conditions 
E: K" —IR+ defined as follows 


E{x) 


d{0 


(1 < P < oo). 


(4.1) 


Let a > 0 and 6 > 1. Throughout this section, we define the function cj) and 
the real number R by (1.7) and (1.8), respectively. It is easy to show that 
R is the unique solution of the equation (f(t) = 1 in the interval [0,r), where 
r = 2/(5 -I-1 -I- a/( 5 — 1)^ -I- 4a). Note that (f is an increasing function which maps 
[0,i?] onto [0,1]. Besides, (j) is quasi-homogeneous of degree 2 on [0,i?]. In the next 
definition, we introduce a sequence of such functions. 


Definition 4.1 We define the sequence (/>jv)?/=o nondecreasing functions 
i/jv: [0, 77] —t [0,1] recursively by setting (j)Q{t) = 1 and 


, __ at^cfNjt) _ 


(4.2) 


where a > 0 and 5 > 1 are constants. 


Proof of the correctness of Definition 4.1 We prove the correctness of the 
definition by induction. For = 0 it is obvious. Assume that for some TV > 0 the 
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function (f>N is well-defined and nondecreasing on [0,i?] and (j)N{R) = 1. We shall 
prove the same for 4>n+i- It follows from the induction hypothesis that 

(1 — t)(l — bt) — > (1 — t)(l — bt) — at^ > Q for all t S [0, i?] (4.3) 


which means that the function (j)]s[+i is well-defined on From (4.2) and the 

induction hypothesis, we deduce that 4>n+i is nondecreasing on [0,R]. From (4.2) 
and (f)pf{R) = 1, we obtain 


4’n+i{R) 


aR?(j)i<i{R) 

(T - .Rj(r-6 R)"{r) 


aR^ 

(1 - R){1 - bR) - aR^ 


cj,{R) = 1. 


This completes the induction and the proof of the correctness of Definition 4.1. □ 


Definition 4.2 For any integer TV > 0, we define the function ip^ : [0, R] —)■ [0, i?] 
as follows 


where the function (j)j^ is defined by Definition 4.1. 


(4.4) 


In the next lemma, we present some properties of the functions and 


Lemma 4.3 Let N > 0. Then: 

(i) (pN is a quasi-homogeneous function of degree 2N on [0,1?]; 

(ii) (fN+iit) < (fit) (fNit) for every t G [0, R]; 

(iii) (fN+iit) < (fNit) for every t G [0,1?]; 

(iv) (/)jv(t) < (fit)^ for every t G [0,1?]; 

(v) pi\f is a gauge function of order 2N -|- 1 on [0,1?]. 


Proof Claim (i) can easily be proved by induction. From (4.2) 


and (4.3), we get 


4>N+lit) 


at^cfiait) 

(1 — t)(l — bt) — at‘^(j)]\[it) 


^ at^(j)Nit) 

~ (1 — t){l — bt) — af^ 


(fit) (f^it) 


which proves (ii). Claim (iii) is a trivial consequence from (ii). Claim (iv) follows 
from (ii) by induction. Claim (v) follows from (i) and the definition ol D 


Lemma 4.4 Let f G IK[z] be a polynomial of degree n > 2, ^ G K" he a root-vector 
of f and N >0. Suppose x G Djq is a vector such that fixi) yf 0 for some i G In- 
(i) If xffT^^\x), then 


f'i^i) _ _ I _ 

/(^*) ^,x,-T^^\x) 


I - CTj 


(4.5) 


where cr^ G K zs defined by 


Ti = ixi-^^)Y^ 




^ {Xi - ^j)ixi - Tj^\x)) 


(4.6) 
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(ii) If X G Djq+i, then 

i Gi 

Proof (i) Taking into account that ^ is a root-vector of /, we get 


(4.7) 


f'i^i) _ 1 


n 1 ^ 

^ (x) 

1 


1 

:37: + E 


Xi - 6 


j¥^i 


-E 


(N), 


1 - (Ti 


Xi {xi - ^j){xi - Tf\x)) Xi 


which proves (4.5). 

(ii) It follows from x € D^+i that 


f\xi) - f{Xi)^ 


1 


^ Xi-Tj^\x) 
Then from (1.13) and (4.5), we obtain 


7^0. 


(4.8) 


T^^~^^\x)-^i = Xi-^i- 


= Xi-^i- 


-E 


fjxi) 

Xi ^i f * \ 

1- {Xi - £.i) 

I - (Ji 1 - cr,: 


which completes the proof. 


□ 


Lemma 4.5 Let f G K[z] be a polynomial of degree n > 2 which has only simple 
zeros in K, ^ S K” he a root-vector of f, N >0 and 1 < p < oo. Suppose x € K" is 
a vector satisfying the following condition 


E{x) 


x-f, 

d{0 


<R = 


2 

6 -I-1 -I- a/( 6 - 1)2 -I- 8a ’ 


(4.9) 


where the function E: K" —)■ IR+ is defined by (4.1), a = (n — and b = 2^/'^. 

Then 


X e Dm, ||TW(a;)-e|| ^ (^7v(if(a;))||x-e|| and E{T^^\x)) < ipN{E{x)). (4.10) 


Proof We shall prove statements by induction on A^. If = 0, then (4.10) holds 
trivially. Assume that (4.10) holds for some N > 0. 

First, we show that x G Dm+i, i.e. x ff T'^^\x) and (4.8) holds for every i S It 
follows from the first inequality in (4.10) that the inequality (3.3) is satisfied with 
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= X and V = T^^\x). Then by Lemma 3.3 and (4.9), we obtain 


n(N), 


1 -. 

x-^ 

1 

66 


\i^-i^>{l-hE{x))d,{^,)>Q (4.11) 


for every j ^ i. Consequently, x # T'^^\x). It remains to prove (4.8) for every i £ In- 
Let i G In he fixed. We shall consider only the non-trivial case f{xi) ^ 0. In this 
case (4.8) is equivalent to 


f(Xi) 


-E 


^0. 


(4.12) 


We define ai by (4.6). It follows from Lemma 4.4(i) that (4.12) is equivalent to 
(Ti 7 ! 1. By Lemma 3.1 with u = x and v = ^ and (4.9), we get 


1 

x-^ 


66 


16-C,I > (i-if(x))d,(6 >0 


(4.13) 


for every j 7 ^ i. From the triangle inequality in K, (4.11), (4.13), induction hypoth¬ 
esis and Holder’s inequality, we get 


W^\ < |a;i-6lE 




< 


< 


1 


Xi - 61 \xi - T^^\x) 


il-E{x))il-bEix)) 6(6 ^ 6(6 


aE{x)ipN{E{x)) 


aE{x)‘^4’N{E{x)) 


(1 - F;(a;))(l - bE{x)) (1 - E{x)){l - bE{x)) ' 
From this, (j)M{E{x)) < 1 and (4.9), we obtain 
aE{x)‘^ 


(4.14) 


Wi\ < 


{l-E{x)){l-bE{x)) 


< 1 


which yields ai ^ 1 and so (4.8) holds. Hence, x G Htv+i- 
Second, we show that the inequalities in (4.10) hold for TV -|- 1. The first inequality 
for -I- 1 is equivalent to 


It}^~''^\x) - 61 < ^n+i(E(x))Ixi - 61 for all i G In 


(4.15) 


Let i G Inhe fixed. If xt = 6 ) then = 6 and so (4.15) becomes an equal¬ 

ity. Suppose Xi 7^ 6- By Lemma 4.4(ii), the triangle inequality in IK and the estimate 
(4.14), we get 


|7;^^+^^(a;)-61 = 11 'E < -I Til 

11 ^i\ 1 \^i\ 

aE{x)'^(j)isi{E{x)) 


< 


(1 — E(x))(l — bE(x)) — aE(xY4>N{E(x)) 

= (j)N+l{E{x))\Xi-^i\ 


\Xi 61 
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which proves (4.15). Dividing both sides of the inequality (4.15) by di{^) and taking 
the p-norm, we obtain 

EiT^^+^\x)) < ^N+iiE{x)) 

which proves that the second inequality in (4.10) holds for + 1. This completes 
the induction and the proof of the lemma. □ 

Now we are ready to state the main result of this section. In the case N = 1 this 
result coincides with Theorem 1.1. 


Theorem 4.6 Let f € IK[z] be a polynomial of degree n > 2 which has only simple 
zeros in K, ^ S K” be a root-vector of f, N > 1 and 1 < p < oo. Suppose x^^^ € K" 
is an initial guess satisfying 


E{x^°^) 




<R = 


P 


2 

5 + 1 + a/(6 - 1)2 + 8a ’ 


(4.16) 


where the function E: K” —>■ IR+ is defined by (4.1), a = (n — and b = 2f^^. 

Then the Ehrlich-type iteration (1.15) is well-defined and converges to f with error 
estimates 


(4.17) 

for all k> 0, where A = (j)N{E(x^^'^)) and the function 4 >n is defined by Defini¬ 
tion 4-L 

Proof We apply Corollary 2.9 to the iteration function : Dn C K" —>■ K" de¬ 
fined by Definition 1.3 and to the function E: K" —>■ K+ defined by (4.1). Let 
J = [0,i?). It follows from the Lemma 4.5, Lemma 4.3(v) and Lemma 2.3 that 
if is a function of initial conditions of with a strict gauge function of order 
r = 2N -I- 1 on J. Since f is a, root-vector of /, then E{^) = 0 G J. It follows from 
Lemma 4.5, that is an iterated contraction at a point f with respect to E 

and with control function (j)^. The fact that x^^'^ is an initial point of follows 
from Lemma 4.5 and Theorem 2.6. Hence, all the assumptions of Corollary 2.9 are 
satisfied, and the statement of Theorem 4.6 follows from it. □ 

Corollary 4.7 Let f G IK[ 2 :] be a polynomial of degree n>2 which has only simple 
zeros inM., f G K" be a root-vector of f, N > 1 and 1 < p < oo. Suppose x^^^ G K" 
is an initial guess satisfying (4.16). Then the Ehrlich-type iteration (1.15) is well- 
defined and converges to f with error estimates 

||^(fc+i) _^|| ^ A^( 2 iv+i)'= ii^(fc) _^|| and ||a;« -^11 ^ -^|| 

(4.18) 


for all k > 0, where A = (j){E{x^^'>)) and (f is a real function defined by (1.7). 
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Proof It follows from Theorem 4.6 and Lemma 4.3(iv). □ 

Let 0 < < 1 be a given number. Solving the equation (j)(t) = h? in the interval 

(0, R), we can reformulate Corollary 4.7 in the following equivalent form. 


Corollary 4.8 Let f € K[z] be a polynomial of degree n > 2 which has n simple 
zeros m K, ^ G IK" be a root-vector of f, N > 1, 1 < p < oo and 0 < h < 1. Suppose 
€ K" is an initial guess which satisfies 


E{x^°'>) 




< Rh 

p 


2 

b + 1 1)2 + Ml + ’ 


(4.19) 


where a= and 5 = 2^/'^. Then the Ehrlich-type method (1.15) is well- 

defined and converges to ^ with error estimates 


||a;(fc+i)-^|| ^ ;, 2 Ar( 2 JV+i)'' ||^(fc)_^|| and ^ (4.20) 


for all fc > 0. 


Remark 4.9 Corollary 4.8 is an improvement of the result of Kjurkchiev and 
Andreev [16] (see Theorem 1.4 above). Suppose that a vector G IK" satisfies 
(1.17). It is easy to show that condition (1.16) is equivalent to the following one: 


0 < c < min 


6 

2(1 + (2n 


_ 25 _ 

l)/i) ’ 3 + 2/i + ^4(3n - 2)/i2 + 4/i + 4n + 1 J ' 


From this, the initial condition (1.17) and 0 < h < 1, we obtain 

-CIloo ^ ch 

S -6 

2h 

3 + 2/i + v/4(3n - 2)/i2 + 4/i + 4n + 1 
2 

3 + v/4(3n - 2) + (4n + l)/h^ 

2 

3 + \/4n — 3 + 4(n — l)/h‘^ 

Therefore, x^^^ satisfies (4.19) with p = oo. Then it follows from Corollary 4.8 
that the Ehrlich-type method (1.15) is well-defined and converges to f with error 
estimates (4.20). From the second estimate in (4.20) and (1.17), we get the estimate 
(1.18) which completes the proof. 


d(0 


5 Local convergence theorem of the second type 

Let / S IK[z] be a polynomial of degree n > 2. We study the convergence of 
the Ehrlich-type method (1.15) with respect to the function of initial conditions 
E: T IR+ defined by 


E(x) 


x-f 

d{x) 


(1 < P < oo). 


(5.1) 
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In the previous section, we introduce the functions (jjN, and the real number 
R with two parameters a > 0 and b> 1. In this section, we consider a special case 
of (j}N, and R when & = 2. In other words, now we define R by 


- 


3 -p ■\/l -p 8(2 


(5.2) 


Furthermore, we define the functions </)n and by Definitions 4.1 and 4.2, respec¬ 
tively, but with 




(1 — t)(l — 2t) — at'^(j)N(t) 


(5.3) 


instead of (4.2), where a > 0 is a constant. 


Definition 5.1 For a given integer iV > 1, we define the increasing function 
Pn- [0,i?] ^ [0,1) by 


^Af(i) 


1 — t — 


(5.4) 


and we define the decreasing function [0,i?] —>■ (0,1] as follows 


1pN{t) = 1 - 2t(l -P ^Ar(t)) 


(1 — t)(l — 2t) — 
1 — t — 


(5.5) 


Proof of the correctness of Definition 5.1 The functions /Sn and are 

well-defined on [0, i?] since 


1 — t — at^(pN-iit) > 1 — t — at^ > 0 for all t £ [0,i?]. 


(5.6) 


The monotonicity of /3n and tpN obvious. It remains to prove that /3jv(i?) < 1 
and ip^iR) > 0. Since (Pn{R) = 1> we obtain 


PNiR) = 


aR^ 


< 1 and ipNiR) = 


(1 - i?)(l - 2R) - aR^ 


l-R-aR^ r.vv / 

which completes the proof of the correctness of Definition 5.1 


> 0 


□ 


Lemma 5.2 Let N > 1. Then: 

(i) (Sn is a quasi-homogeneous of degree 2N on [0, i?]; 

(ii) /3Ar(t) = (pN{t)%pN{t) for every t £ [0,i?]; 

(hi) ^Ar+i(t) < /3Ar(f) for every t £ [0,i?]; 

(iv) '4>N+i{t) > ipN{t) for every t £ [0,i?]. 

Proof The function can be presented in the form piN{t) = t^(pM-i{t)^{t)^ where 
$(t) = a/{1 — t — at^(pM-i{t)). Therefore, /3jv is quasi-homogeneous of degree 2N 
on [0, i?] since it is a product of three quasi-homogeneous functions on [0, R] of 
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degree 2, 2N — 2 and 0. From the definitions of the functions ipjsf, and /Sn, we 
get 


(pN{t)tpN{t) = 


(1 — t)(l — 2t) — 

(1 — <)(1 — 2t) — 1 — f — at^4>N-i(t) 




Claim (iii) follows from Lemma 4.3(iii) and (5.4). Claim (iv) follows from (iii) and 
(5.5). □ 


Lemma 5.3 Let f S K[ 2 ;] be a polynomial of degree n> 2 which has only simple 
zeros in K, ^ G K" a root-vector of f, N > 1 and 1 < p < oo. Suppose x € K" is a 
vector with distinct components such that 


E{x) 


d{x) 


<R = 


2 


(5.7) 


where the function if: D —>■ K_|_ is defined by (5.1) and a = (n— 1)^/'^. Then f has 
only simple zeros in K, 


cr G Dn, ||TW(cr)-C|| ^ ^^(F;(cc))||cr-C|| and E{T<^^\x)) < ^n{E{x)). (5.8) 


Besides, the vector T^^^x) has pairwise distinct components. 


Proof It follows from (5.7) and R < 1/2 that E{x) < 1/2. Then it follows from 
Lemma 3.2 that the vector ^ has distinct components, which means that / has only 
simple zeros in K. We divide the proof into two steps. 


Step 1 In this step, we prove x G Dn and the first inequality in (5.8) by induction 
on If = 1, the proof of the claims can be found in [11]. Assume that x & Dm 
and the first inequality in (5.8) hold for some A^ > 1. 

First we show that x G Dm+i i-G. xfpT^^\x) and (4.8) holds for every i G It 
follows from the first inequality in (5.8) that (3.5) holds with v = x, u = T^^\x) 
and a = 1. Therefore by Lemma 3.4, (5.7) and R < 1/2, we obtain 


x,-tI^\x)\ > 1-2 


x-^ 


d{x) 


\xi — Xj\ > (1 — 2E{x)) dj{x) > 0 (5.9) 


for every j yf i. Consequently, x ff {x). It remains to prove (4.8) for every i G In- 
Let i G Ira be fixed. We shall consider only the non-trivial case f{xi) 0. In this 
case (4.8) is equivalent to (4.12). On the other hand, it follows from Lemma 4.4(i) 
that (4.12) is equivalent to Cj yf 1, where ai is defined by (4.6). By Lemma 3.1 with 
u = ^ and v = X and (5.7), we get 


1 - 

x-f. 

1 

d{x) 


\Xi fj I ^ 


P 


\xi-Xj\ = {1-E{x))\xi-Xj\ > {1-E{x)) di{x) > 0, 

(5.10) 
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for every j ^ i. Hence, we obtain From induction hypothesis, we get 




(5.11) 


Combining the triangle inequality in K, (5.10), (5.9) and (5.11), we obtain 


kil < 


< 


< 


^ I ^ I I M 

j^^ \xi - ^j\\xi - Tj '(a;) I 

{I-E{x)){l-2E{x)) di{x) ^ dj{x) 

Pn{E{x)) \xj - ^j\ 

{I - E{x)){l - 2E{x)) di{x) ^ dj{x) 


which, using Holder’s inequality, yields 

^ aE{xY I3 m{E{x)) 
{l-E{xm-2E{x))- 

From this and (5.7), we deduce 

. . ^ aE{x)^ 

- (l-F;(a:))(l-2F;(a:)) ^ 


(5.12) 


which yields ai ^ 1, and so (4.12) holds. Thus we prove that x G Hat+i. 

Now we have to prove that the first inequality in (5.8) is satisfied for + 1, which 
is equivalent to 


\Tj;^~''^\x) - ^i\ < PN+i{E{x))\xi - ^i| for all i G (5.13) 

Let i G In be fixed. If Xi = ^i, then (x) = and the inequality (5.13) becomes 

an equality. Suppose Xi ^ It follows from Lemma 4.4(ii), the triangle inequality 
in K and the estimate (5.12) that 




(a;)-61 = 


< 


k»l 

|1 - cr*| 


\Xi 


61 


< 


1 - \<7i\ 


\x^ 


61 


aE{xY I3 n{E{x)) 

(1 — E{x)){l — 2E{x)) — aE{xYldpf{E(x)) 


61- 


From this inequality. Lemma 5.2(ii), YNit) 6 !> (5.3) and Lemma 5.2(iv), we obtain 


< 

< 


_ aE{xY(j)N{E{x))YN{E{x)) _ _ 

(1 - E{x)){l - 2E{x)) - aE{xY(j)N{E[x)) 

(()7V+i(^^(a:))V’Af(^^(a:))|a;i - 61 

(('7V+i(^^(a:))'!/’Af+i(^^(a:))|a;j - 61 = l3N+i{E{x))\xi - 6| 


which proves (5.13). This completes the induction. 
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step 2 In this step we prove the second inequality in (5.8) and that has 

distinct components. First inequality in (5.8) allow us to apply Lemma 3.5 with 
u = V = X and a = Pn{E{x)). By Lemma 3.5 and (5.5), we deduce 

\T^^\x) -T^^\x)\ > {l-2^/'^E{x){l + PNiE{x))))\x^-Xj\ > ^N{E{x))\xi-Xj\. 

By taking the minimum over all j G In such that j ^ i, we obtain 

d,(TW(cr)) > ijN{E{x))d,{x) > 0 (5.14) 

which implies that T^^\x) has distinct components. It follows from (5.11), (5.14) 
and Lemma 5.2(ii) that 

\Ti^\x) -^i\ I3n{E{x)) |x, - ^i\ ^ -C»l 

di{T(^'){x)) ~ tj}N{E{x)) di{x) ^ di{x) 

By taking the p-norm, we obtain 

E{T^^\x)) < cj>N{E{x)) E{x) = vn{E{x)) 

which proves the second inequality in (5.8). This completes the proof. □ 

Now we are able to state the main result of this section. In the case when = 1 

and p = oo this result reduces to Theorem 1.2. 


Theorem 5.4 Let f G IK[z] be a polynomial of degree n > 2 which splits over K, 
^ G K” be a root-vector of f, N > 1 and 1 < p < oo. Suppose x^^^ G K" is an initial 
guess with distinct components such that 


E{x^°'>) 


d(a:(°)) 


< R = 


3 +Vl + 8a’ 


(5.15) 


where the function E is defined by (5.1) and a = {n — 1)^/"^. Then f has only simple 
zeros in K and the Ehrlich-type iteration (1.15) is well-defined and converges to f 
with error estimates 


(5.16) 

for all k> 0, where A = (j)j^{E{x^'^'^)), 9 = ipt^{E{x^^'^)). Moreover, the method is 
convergent with order 2N + 1 provided that E{x^^'^) < R. 

Proof We apply Theorem 2.8 to the iteration function : Dn C K" —>■ K” to¬ 
gether with the function E: Dj^ —>■ K_|_ defined by (5.1). 

It follows from Lemma 5.3 and Lemma 4.3(v) that if is a function of initial 
conditions of with gauge function pjv of order r = 2N -|- 1 on the interval 

J= [0,ii]. 
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From Lemma 5.3, we get that is an iterated contraction at ^ with respect to 
E and with control function /?jv. Also, it is easy to see that the functions /3pf, 

'tpN and ifiN have the properties (2.5), (2.6) and (2.7). 

It follows from Lemma 5.3 that S D^- According to Theorem 2.6 to prove 
that is an initial point of it is sufficient to prove that 

x€Dn and E{x)€J => T^^^Hx) € Dn ■ (5.17) 

From X € D^, we have T^^\x) G K". By Lemma 5.3, has distinct compo- 
nents and E{T^^\x)) < ipn{E{x)). The last inequality yields E{T^^\x)) G J since 
ip: J ^ J and E{x) G J. Thus we have both T^^\x) G D and (x)) G J. Ap¬ 

plying Lemma 5.3 to the vector T^^\x), we get T^^\x) G Dj^ which proves (5.17). 
Therefore, is an initial point of 

Now the statement of Theorem 5.4 follows from Theorem 2.8. □ 


6 Semilocal convergence theorem 

In this section we establish semilocal convergence theorems for Ehrlich-type meth¬ 
ods (1.15) for finding all zeros of a polynomial simultaneously. We study the conver¬ 
gence of these methods with respect to the function of initial conditions E: D —>■ M+ 
defined by 


Ef{x) 


Wfjx) 

d{x) 


(1 < P < c»). 


( 6 . 1 ) 


Recently Proinov [22] has shown that there is a relationship between local and 
semilocal theorems for simultaneous root-finding methods. It turns out that from 
any local convergence theorem for a simultaneous method one can obtain as a conse¬ 
quence a semilocal theorem for the same method. In particular, from Theorem 4.6 
we can obtain a semilocal convergence theorem for Ehrlich-type methods (1.15) 
under computationally verifiable initial conditions. For this purpose we need the 
following result. 


Theorem 6.1 (Proinov [22]) Let / G K[z] he a polynomial of degree n> 2. Sup¬ 
pose X G K" is an initial guess with distinct components such that 


Wfjx) 

d{x) 


R{1 - R) 

1 -1- (a - l)i? 


( 6 . 2 ) 


for some 1 < p < oo and 0 < R < 1/(1 -l- \/a), where a = {n — 1)^/*. In the case 
n = 2 and p = oo we assume that inequality in (6.2) is strict. Then f has only 
simple zeros in IK and there exists a root-vector ^ G K" of f such that 


||x-^|| ^ a(F;/(a:)) ||lF(a;)|| and 


d{x) 


< R, 

p 


(6.3) 


where the real function a is defined by 


a{f) 


2 

1 — (a — l)t -I- a/(1 — (a — f)tY — 4t 


(6.4) 
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If the inequality (6.2) is strict, then the second inequality in (6.3) is strict too. 


Now, we are ready to state and prove the main result of this paper. 


Theorem 6.2 Let f G be a polynomial of degree n > 2, N > 1, 1 < p < oo. 
Suppose G K" is an initial guess with distinct components such that 


Efix^^^) 


lT(a;(°)) 

(i(a;(°)) 


8 

p ^ (3 + v/rT^)2 ’ 


(6.5) 


where the function Ef is defined by (6.1) and a = {n — 1)^/'^. Then f has only simple 
zeros in K and the Ehrlich-type iteration (1.15) is well-defined and converges to a 
root-vector f of f with order of convergence 2N + 1 and with a posteriori error 
estimate 


||ir«-CHa(i?/(^('=^))l|W^/(ir«)||, (6.6) 

for all k > 0 such that Ejfx^^'^) < 8/(3 + Vl + Safi, where the function a is defined 
by (6.4). 


Proof Let us define R by (5.2). It is easy to calculate that R < 1/(1 + ^/a) and 


i?(l - R) 2(1 + Vl + 8a) 8 

1 + (a - l)i? “ (3 + v/rT8a)(l + 2a + VI + 8a) “ (3 + v/TTSo)^ ' 


Therefore, (6.5) can be written in the form 


lT(a;W) 

d(xi°)) 


R{1 - R) 

1 + (a - 1)R ■ 


Then it follows from Theorem 6.1 that / has only simple zeros in K and there exists 
a root-vector f G K" of / such that 


d{x^^'>) 


< R. 


Now Theorem 5.4 implies that the Ehrlich-type iteration (1.15) converges to ^ with 
order of convergence 2N -P 1. It remains to prove the error estimate (6.6). Suppose 
that for some fc > 0, 


VE/(a;('=)) 

d{x^’^'>) 


R{l-R) 
l + ia-l)R ■ 


(6.7) 


Then it follows from Theorem 6.1 that there exists a root-vector ij G K” of / such 
that 


— ryll ^ a{Ef{x^'^'^)) ||lE/(a;^^^)|| and 


2 ;(fe) _ jy 

d(a;('')) 


< R. 

p 


( 6 . 8 ) 
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From the second inequality in (6.8) and Theorem 5.4, we conclude that the Ehrlich- 
type iteration (1.15) converges to rj. By the uniqueness of the limit, we get rj = ^. 
Therefore, the error estimate (6.6) follows from the first inequality in (6.8). This 
completes the proof. □ 


Setting p = oo in Theorem 6.2, we obtain the following result. 


Corollary 6.3 Let f G K[z] be a polynomial of degree n >2 and N > 1. Suppose 
S IK is an initial guess with distinct components such that 


d(a;(°)) 


8 

(3 + v/8h:^j2' 


(6.9) 


Then f has only simple zeros in IK and the Ehrlich-type iteration (1.15) is well- 
defined and converges to a root-vector ^ of f with order of convergence 2N + 1 and 
with error estimate (6.6) for p = oo. 


Setting p = 1 in Theorem 6.2 we obtain the following result. 


Corollary 6.4 Let f G IK[z] be a polynomial of degree n > 2 and N > 1. Suppose 
G IK" is an initial guess with distinct components such that 


lT(a;W) 

d(a;(°)) 



( 6 . 10 ) 


Then f has only simple zeros in IK and the Ehrlich-type iteration (1.15) is well- 
defined and converges with order 2N -\- 1 to a root-vector f of f with error estimate 
(6.6) for p = 1. 


7 Numerical examples 

In this section, we present several numerical examples to show some applications of 
Theorem 6.2. Let / G C[z] be a polynomial of degree n>2 and let x^^'^ G C" be an 
initial guess. We show that Theorem 6.2 can be used: 

• to prove numerically that / has only simple zeros; 

• to prove numerically that iVth Ehrlich-type iteration (1.15) starting from x^^'^ 
is well-defined and converges with order 2N -|- 1 to a root-vector of /; 

• to guarantee the desired accuracy when calculating the roots of / via A^th 
Ehrlich-type method. 

In the examples below, we use the function of initial conditions Ef. T) ^ K+ 
defined by 


Ef{x) 


Wfjx) 

d{x) 


(7.1) 


where Wf is the Weierstrass correction defined by (1.3). We consider only the case 
p = oo since the other cases are similar. 
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Also, we use the real function a defined by 

2 


i{t) = 


1 — (n — 2)t + a/(1 — (n — 2 )i )2 — 

It follows from Theorem 6.2 that if there exists an integer m > 0 such that 


(7.2) 


(3 + y/8n - 7)2 ’ 


(7.3) 


then / has only simple zeros and the Ehrlich-type iteration (1.15) is well-defined 
and converges to a root-vector ^ of / with order of convergence 2N + 1. Besides, 
for all fc > TO such that 


E/(a;('=)) <^ = 


(3 -f V 8 n - 7)2 
the following a posteriori error estimate holds: 

- ^lloo < Efc, where et = a{Ef{x^'"'>)) \\Wf{x^'"^)\\ao ■ 


(7.4) 


(7.5) 


In the examples, we apply the Ehrlich-type methods (1.15) for some N > 1 using 
the following stopping criterion: 


Efc < 10 and Ef{x^^'^) < {k >m). (7.6) 

For given N we calculate the smallest to > 0 which satisfies the convergence 
condition (7.3), the smallest A: > to for which the stopping criterion (7.6) is satisfied, 
as well as the value of Sk for the last k. 

In Table 2 the values of iterations are given to 15 decimal places. The values of 
other quantities Ef{x^'^'>), etc.) are given to 6 decimal places. 

Example 7.1 We consider the polynomial 


f{z) = z^-l 


and the initial guess 

= (0.5 -f 0.5f, -1.36 -f 0.42f, -0.25 -f 1.28q 0.46 - 1.37f) 

which are taken from Zhang et al. [23]. We have = 0.125 and E{x^^'^) = 0.506619. 
The results for this example are presented in Table 1. For example, we can see that 
for = 10 at the first iteration we have proved that the Ehrlich-type method 
converges with order of convergence 21 and that at the second iteration we have 
calculated the zeros / with accuracy less than 10“^^^. Moreover, at the next itera¬ 
tion we obtain the zeros of / with accuracy less than 10“^®®^. Also, we can see that 
for N = 100 at the second iteration we have obtained the zeros of / with accuracy 
less than 10 “ 
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Table 1 Values of m, k and ek for Example 7.1 {JX = 0.125) 


N 

m 




k 



Sk + 1 

1 

2 

0.010032 

1.457548 X 

10-2 

4 

4.385760 

X 

10-21 

8.919073 

X 

10-62 

2 

1 

0.067725 

1.242914 X 

10-1 

3 

1.347060 

X 

10-28 

7.284576 

X 

10-192 

3 

1 

0.015716 

2.300541 X 

10-2 

3 

1.825502 

X 

10-106 

5.054741 

X 

10-744 

4 

1 

0.002730 

3.887455 x 

10-2 

2 

1.330837 

X 

10-25 

3.543773 

X 

10-220 

5 

1 

0.001215 

1.722883 X 

10-2 

2 

4.720064 

X 

10-27 

2.999643 

X 

10-407 

6 

1 

0.000206 

2.927439 x 

10-4 

2 

1.060096 

X 

10-50 

5.523501 

X 

10-657 

7 

1 

0.000081 

1.155284 X 

10-4 

2 

6.261239 

X 

10-67 

3.252761 

X 

;j^Q-1002 

8 

1 

0.000014 

1.986052 X 

10-2 

2 

6.080606 

X 

10-85 

3.570038 

X 

10-1439 

9 

1 

0.000005 

7.910775 X 

10-2 

2 

1.309022 

X 

10-105 

1.170454 

X 

]^q-2002 

10 

1 

0.000000 

1.366899 X 

10-2 

2 

4.301615 

X 

10-128 

8.477451 

X 

^q-2683 

100 

1 

0.000000 

1.820743 X 

10-57 

1 

1.820743 

X 

10-57 

3.460397 

X 

10-11451 


In Table 2, we present numerical results for Example 7.1 in the case N = 10. 


Table 2 Numerical results for Example 7.1 in the case JV = 10 


rh li. ^ 

0 (15 +051 

1 1.000000380419496 + 0.000000816235730i 

2 1.000000000000000 + O.OOOOOOOOOOOOOOOf 


fu tJL ^ 

0 -0.25 + 1.28i 

1 0.000000277962637 + 0.999999578393062i 

2 0.000000000000000 + l.OOOOOOOOOOOOOOOi 


+ 0.421 

-1.000000220051461 - 0.000000495915480f 
-1.000000000000000 + O.OOOOOOOOOOOOOOOf 


_ 1^4 _ 

0.46 - 1.37* 

-0.000000314533436 - 0.999998669784542i 
0.000000000000000 - l.OOOOOOOOOOOOOOOi 


Example 7.2 We consider the polynomial 


f(z) = Z^^+Z^^ + l 


and Aberth’s initial approximation € C" given by (see Aberth [8] and Petkovic 
et al. [24]): 

=-—+roexp(i0^), = iy = l,...,n, (7.7) 

n n \ 2 J 

where ai = 1, rp = 2 and n = 15. We have ^ = 0.043061 and E{x^^'^) — 0.179999. 
The results for this example are presented in Table 3. For example, we can see that 
for = 30 at the third iteration we have obtained the zeros of / with accuracy less 
than 10“^“^®. Moreover, at the next iteration we get the zeros of / with accuracy 
less than 10“^®^°®. 


Table 3 Values of m, k and for Example 7.2 = 0.043061) 


N 

m 




k 

£k 


£fc+l 

1 

6 

0.036897 

3.187918 X 

10-2 

9 

3.967908 

X 

10-36 

5.304009 

X 

10-106 

2 

5 

0.000003 

1.182714 X 

10-6 

6 

6.112531 

X 

10-28 

2.230412 

X 

10-124 

3 

4 

0.000064 

2.475020 X 

10-5 

5 

2.446120 

X 

10-29 

2.722168 

X 

10-197 

4 

4 

0.000000 

1.550670 X 

10-11 

5 

3.838741 

X 

10-92 

1.589981 

X 

10-827 

5 

3 

0.005793 

2.415745 X 

10-2 

4 

9.532339 

X 

10-24 

8.487351 

X 

10-248 

6 

3 

0.000293 

1.127450 X 

10-4 

4 

9.565008 

X 

10-45 

1.725858 

X 

10-565 

7 

3 

0.000005 

2.173198 X 

10-6 

4 

4.018844 

X 

10-77 

6.737932 

X 

10-1128 

8 

3 

0.000000 

1.562375 X 

10-2 

4 

1.162424 

X 

10-122 

1.291370 

X 

]^q-2080 

9 

3 

0.000000 

4.092421 X 

10-11 

4 

4.245137 

X 

10-187 

1.373908 

X 

Iq-3530 

10 

3 

0.000000 

3.904607 X 

10-14 

4 

4.643262 

X 

10-270 

2.543247 

X 

]^q-5644 

30 

2 

0.000055 

2.129417 X 

10-5 

3 

5.721566 

X 

10-249 

2.377023 

X 

j^q-15106 















Proinov and Vasileva 


Page 24 of 26 


Example 7.3 We consider the Wilkinson polynomial ([25]) 


20 


f{z) = Y[{z-j)= - 120^1^ + ... + 2 432 902 008 176 640 000. 


and Abert’s initial approximation (7.7) with oi = —120, tq = 20 and n = 20. We 
have ^ = 0.033867 and = 0.344409. The results foe Example 7.3 are shown 

in Table 4. For example, we for N = 100 at the seventh iteration we get the zeros 
of / with accuracy less than 10“^^^^®. 


Table 4 Values of m, k and for Example 7.3 = 0.033867 ) 


N 

m 




k 



£fc+l 

1 

18 

0.000060 

6.095859 X 

10-5 

20 

1.620028 

X 

10-38 

4.276235 

X 

10-114 

2 

12 

0.015335 

2.153155 X 

10-2 

14 

1.095084 

X 

10-46 

1.779476 

X 

10-250 

3 

10 

0.018005 

2.769333 x 

10-2 

12 

8.917532 

X 

10-86 

4.482714 

X 

10-596 

4 

9 

0.005514 

6.130790 X 

10-5 

10 

4.221856 

X 

10-21 

7.250879 

X 

10-184 

5 

9 

0.000000 

1.159694 X 

10-55 

10 

5.021359 

X 

10-165 

5.118016 

X 

j^q-1808 

6 

8 

0.000237 

2.386016 X 

10-* 

9 

8.455240 

X 

10-48 

1.280870 

X 

10-612 

7 

8 

0.000000 

2.723047 X 

10-17 

8 

2.723047 

X 

10-17 

8.926059 

X 

10-249 

8 

7 

0.018995 

2.934241 X 

10-2 

8 

2.885374 

X 

10-30 

4.152134 

X 

10-505 

9 

7 

0.002180 

2.274734 x 

10-5 

8 

3.792876 

X 

10-51 

1.140751 

X 

10-958 

10 

7 

0.000000 

5.185525 X 

10-7 

8 

1.620086 

X 

10-152 

2.936276 

X 

]^q-2768 

30 

5 

0.000181 

1.821419 X 

10-^ 

6 

1.395923 

X 

10-226 

1.902920 

X 

10-15777 


In the Figure 1, we present the trajectories of approximations generated by the 
method (1.15) for N = 30 after 6 iterations. 



Figure 1 Trajectories of approximations for the Wilkinson polynomial f{z) = — j) 
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Example 7.4 We consider the polynomial 


f{z) = - 1. 


In this example we nse Abert’s initial approximation (7.7) with oi = 0, rg = 2 and 
n = 40. We have M = 0.018685, = 0.159318. The resnlts for Example 7.4 

can be seen in Table 5. 


Table 5 Values of m, k and ej, for Example 7.4 = 0.018685) 


N 

m 




k 



£)c+l 

1 

15 

0.007235 

1.588799 X 

10-3 

17 

1.057241 

X 

10-18 

1.574672 

X 

10-52 

2 

11 

0.000001 

1.731641 X 

10-" 

12 

2.763909 

X 

10-30 

2.863869 

X 

10-144 

3 

9 

0.000026 

4.171842 X 

10-® 

10 

5.167701 

X 

10-32 

2.328540 

X 

10-213 

4 

8 

0.000032 

5.141616 X 

10-6 

9 

7.830010 

X 

10-40 

3.487627 

X 

10-344 

5 

7 

0.010766 

2.954474 X 

10-3 

8 

1.468181 

X 

10-20 

2.870206 

X 

10-208 

6 

7 

0.000002 

4.201055 X 

10-" 

8 

7.096655 

X 

10-71 

6.481892 

X 

10-®oo 

7 

7 

0.000000 

9.445503 X 

10-3® 

8 

3.169914 

X 

10-196 

2.445585 

X 

]^q-2918 

8 

6 

0.010675 

2.911647 X 

10-3 

7 

8.218559 

X 

10-31 

3.538870 

X 

10-495 

9 

6 

0.000281 

4.462548 X 

10-® 

7 

2.324176 

X 

10-64 

1.205364 

X 

10-1190 

10 

6 

0.000000 

1.231259 X 

10-7 

7 

1.392265 

X 

10-124 

1.840079 

X 

]^q-2580 

30 

5 

0.000000 

2.416285 X 

10-34 

5 

2.416285 

X 

10-34 

1.294365 

X 

10-1987 


In the Figure 2, we present the trajectories of approximations generated by the 
method (1.15) for N = 30 after 5 iterations. 


iv = 30 = 5) 



Figure 2 Trajectories of approximations for the polynomial f{z) = — 1 
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